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Abstract 
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We present a fully automated method for top-down segmentation of the pul- 
^ monary arterial tree in low-dose thoracic CT images. The main basal pul- 

monary arteries are identified near the lung hilum by searching for candidate 
^ vessels adjacent to known airways, identified by our previously reported air- 

, O, way segmentation method. Model cylinders are iteratively fit to the vessels 

to track them into the lungs. Vessel bifurcations are detected by measuring 
^ the rate of change of vessel radii, and child vessels are segmented by initiating 

O new trackers at bifurcation points. Validation is accomplished using our novel 

^ sparse surface (SS) evaluation metric. The SS metric was designed to quan- 

V/^ tify the magnitude of the segmentation error per vessel while significantly 

\Q decreasing the manual marking burden for the human user. A total of 210 

O arteries and 205 veins were manually marked across seven test cases. 134/210 

^ arteries were correctly segmented, with a specificity for arteries of 90%, and 

average segmentation error of 0.15 mm. This fully-automated segmentation 
is a promising method for improving lung nodule detection in low-dose CT 
screening scans, by separating vessels from surrounding iso-intensity objects. 
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1. Introduction 

The pulmonary arteries comprise the vessels that supply deoxygenated 
blood to the lungs for gas exchange. Computer segmentation of the pul- 
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monary vessels in computed tomography (CT) images has been proposed as 
a method for improving automated lung nodule detection in low-dose lung 
cancer screening scans [Ij. Additionally, computer segmentation of the pul- 
monary vessels has been explored as a method to reduce false positives in 
pulmonary embolism detection by radiologists [2j, and improve pulmonary 
artery size measurements for pulmonary arterial hypertension patients [3]. 

The pulmonary vasculature is a highly complex branching structure that 
presents significant challenges towards robust segmentation. Normal anatom- 
ical variability, as well as anatomic abnormality from disease, limits the 
amount of a priori information that can be used to develop a robust seg- 
mentation algorithm. In low-dose CT lung cancer screening scans, there is 
the added challenge of a high signal-to-noise ratio compared with standard 
diagnostic CT scans. Additionally, lung cancer screening scans do not use 
intravenous radiocontrast, providing low levels of contrast between the dif- 
ferent types of dense (non-fat) soft tissue, such as arteries, veins, and lung 
nodules. 

We present a model-based automated method for segmenting the pul- 
monary arteries in low-dose CT scans and separating them from surrounding 
isointensity objects. Seed-points for initiating a cylindrical vessel tracker are 
automatically located in the basal pulmonary arteries by searching near air- 
ways identified using our automated airway segmentation algorithm [4J . The 
vessel tracker iteratively follows the arteries into the lung parenchyma, and 
vessel branching is detected by means of a bifurcation detector based on the 
rate of change of the vessel radii. 

Attempts to segment the pulmonary vasculature in CT have included 
region-growing [6l [7] , vessel enhancement with a Hessian matrix [SJ O [lO] , 
and fuzzy connectivity [llj. Similar techniques have also been used to seg- 
ment arteries and veins separately. Lei et al. [12] used multiseeded fuzzy 
connectivity to separate a confiated A/V tree in contrast-enhanced mag- 
netic resonance angiography of the abdomen. Bemmel et al implemented 
a level-set algorithm to separate arteries and veins in magnetic-resonance 
images ^J. Recently, Saha et al. presented an approach for separating 
fused isointensity objects, including pulmonary arteries and veins, which uti- 
lized multiscale morphologic filtering, and fuzzy connectivity maps generated 
from manual seed points. 

Although previously-proposed pulmonary artery segmentation methods 
require the use of manually selected seed-points, in this work the segmenta- 
tion is completely automated. This serves to eliminate inter-user differences. 
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and would allow the algorithm to be included in fully-automated lung nodule 
detection systems. Additionally, the method presented here utilizes a cylin- 
drical vessel tracker, which reconstructs vessels by progressively fitting model 
cylinders. Using a model-fitting approach on low-dose scans is advantageous 
because it can identify differences in geometry between two isointensity struc- 
tures, even when local amounts of noise are high. 

We also present a novel evaluation metric, the sparse surface (SS) metric, 
for use with segmentations of vascular trees. The SS metric was designed 
to quantify both the overall extent (rate of correctly segmented vessels) and 
accuracy per identified vessel in the overall segmentation. The validation 
scheme for the SS metric requires a manageably sparse number of manual 
markings, increasing the number of cases that experimental algorithms can 
be validated on in a limited amount of time. 

2. Materials and methods 

2.1. Preliminary steps 

Three preliminary steps were performed to reduce noise and decrease the 
signal from unwanted structures. A 3 x 3 median filter was first applied to 
reduce the non-Gaussian noise inherent in low-dose scans. The scans were 
then thresholded at -400 HU to separate the soft tissue (0 HU) from the lung 
parenchyma (-900 HU) within the lungs. Structures such as bone (200 HU) 
and fat (50 HU) were not eliminated after thresholding. However, because 
these structures do not appear in the lungs, they are highly unlikely to aflFect 
the results of the segmentation. 

We used the results of our previously reported airway segmentation [15] to 
remove the airway walls, a potential source of error in vessel segmentation. 
The airway was enlarged by dilation with a 2 mm isotropic kernel. The 
dilated airway was then subtracted from the CT image to remove the airway 
walls. Figure [T] illustrates the application of this operation to a 3rd generation 
airway/artery segment. 

2.2. Automated seed point detection 

The airways and pulmonary arteries run parallel and adjacent to each 
other in units known as bronchopulmonary segments. This symmetry pro- 
vides the key for our automated pulmonary artery seed point detector. The 
airway lumen is automatically segmented and labelled using a region-growing 
method with an automatically generated seed point in the trachea [15j. Using 
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Original Image Median Filtered Thresholded Airway Walls 

Removed 

Figure 1: Removal of the airway walls by dilation of the airway segmentation and sub- 
traction from the thresholded CT image. 

the segmented airway, a curved cylindrical region of interest is constructed 
around individual airway branches to limit the search space for the pul- 
monary arteries. A cubic spline S{n) is fit to the airway, and the scan is 
resampled to construct cutting planes orthogonal to the spline, spaced at 
constant intervals. Figure ^ dia grams the processes of generating the 3D 
ROI, and shows a single frame for illustration. 




Initial opening: Final opening: 



Artery not separated Artery separated 

Figure 2: a) Diagram of the airway ROI construction process and a single frame from the 
output, b) Progressive morphological opening operation to isolate the pulmonary artery, 
c) Complete process diagram for the identification of seed points within the pulmonary 
artery. 

The algorithm for locating the pulmonary arteries within the reconstructed 
ROI is outlined in Figure The ROI is thresholded at -400 HU, and the 
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soft tissue structures are separated using a 3D morphological opening oper- 
ation to separate the heart from the pulmonary arteries (Figure [2}3). The 
resulting objects are labeled with a 3D connected component algorithm. 

Each unique object is then checked to see if it is a possible candidate for 
the pulmonary artery. Objects which extend past 4 cm from the center of the 
airway contain areas of the heart and are eliminated. Of the remaining ob- 
jects after the elimination, the component that is closest to the airway center 
is identified as the pulmonary artery, as predicted by the bronchopulmonary 
model. If no components are found, the morphological operation opening 
was too small to isolate the artery. In this case, the size of the kernel is in- 
creased by one voxel and the process is repeated until the pulmonary artery 
is identified. 

The location of the isolated pulmonary artery is then defined as the cen- 
troid of the artery, in the ROI frame most proximal to the heart. This ensures 
that vessel tracking begins as close to the start of the vessel as possible, in- 
creasing the extent of the segmentation. The orientation of the artery was 
found by computing the centroid of the artery in the most distal frame it 
appears in the ROI, and taking the vector from the proximal centroid to the 
distal centroid as the initial artery orientation. 

2.3. Vessel tracking 

A single vessel can be modeled as a series of discrete cylinders based 
on the observation that they are approximately cylindrical in shape. Vessel 
tracking is the task of recreating the cylindrical vessel by progressively fitting 
model cylinders in the CT space. 

The similarity metric between a model cylinder and the CT scan is con- 
structed from the cylinder C and set V of voxels in the thresholded image. 
The total similarity score S is the sum of scores Si of all voxels Vi in C: 



The strict penalty of —5 was selected ad hoc to prevent the cylinder tracker 
from jumping between adjacent soft tissue structures model cylinders that 
span areas of lung parenchyma are highly penalized. The similarity metric 
is not normalized to the size of the cylinder, which rewards cylinders that 
capture more of the desired vessel. 




where 




(1) 



5 



The iterative vessel tracking process is outlined in Figure [3^. Beginning 
at seed point Xq and direction do, model cylinders of varying orientations and 
radii are fit to the vessel. The initial radius Tq is optimized on the interval 
[2 mm, 12 mm], to cover both abnormally large and small vessels. The 
orientation search space is limited to the solid angle = 27r, which allows 
for sharp turns in vessels. To limit the size of the model parameter space, 
the cylinder height is set to a fixed value of ho. The segmented vessel is then 
the set of soft-tissue pixels contained within the best fitting model cylinder. 

After the initial model cylinder is fit, the process repeats for a new initial 
position xi, shown in Figure (3)3. For the t^^ iteration of cylinder fitting, the 
initial point Xt is defined as: 

Xt = (Astep)rft-l + Xt-i (2) 

where rft-i is the normalized direction of the previous model cylinder, and 
Astep is the distance the next cylinder is progressed from the previous one. 
Because the radius is expected to vary slowly between model cylinders, the 
radius search space is limited on: 

r, G [0.7r,_i,1.3r,_i] (3) 



Seed point xo ^ 
and direction do 




End vessel 




Find and 




track children 



Figure 3: a) Process diagram for iterative vessel tracking. Starting from an initial seed 
point xq and direction (io, cylinders are iteratively fit to the image, recreating the vessel. 
New vessels are initiated at points of bifurcation, b) Diagram showing the progressive 
fitting of model cylinders. Several model cylinders of varying radii and orientation are 
checked for their similarity to the vessel. 
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Iterative vessel tracking terminates when the tracker either fails to find a 
strong match or a bifurcation is detected. A strong match is defined as having 
at least 50% soft-tissue voxels within the model cylinder. This criterion fails 
to be met when either the vessel becomes too small to give a strong signal, 
as in the case of higher generational subsegmental arteries, or the tracker has 
fallen off the vessel. 

2.4' Bifurcation detection 

To initiate tracking of child vessels, a method was developed to search for 
child branches and initiate new vessel trackers. Automated bifurcation de- 
tection is performed by comparing the actual vessel geometry to a bifurcation 
model based on the following observations: 

1. The radii of the child vessels is smaller than the radii of their parent 
vessel. 

2. The angle between child vessels is no less than 30^ and no more than 
90^ 

3. The child vessels curve only slowly away from the parent vessel. 

As the parent tracker encounters a bifurcation, the next model cylinder 
will fit onto the child vessel with a smaller radius as shown in Figure [4)3. 
Using this model, a bifurcation candidate is generated when the ratio of the 
current radius to the radius at I ho upstream, ro/r^o, falls below the radius 

change threshold ^radius- 

The process for finding child vessels at bifurcation candidate points is 
outlined in Figure [3^. At candidate point x^, a new child vessel (child 1) is 
initiated at Xt in the direction dt. The next step is to search for the second 
child vessel. Observations 2 and 3 of the bifurcation model indicate where to 
look for bifurcation candidates. Let the orientation of the parent vessel at the 
bifurcation be rfp, and the orientation of children 1 and 2 at the bifurcation 
be di and ^2 respectively. Child vessel 2 will satisfy the bifurcation model if: 

cos ^ < (^di • d2^ < cos ^ and • dp^ < cos ^ (4) 

The search space that satisfies these conditions is illustrated in Figure [4} 

In the case of a false positive in the bifurcation detection, child 1 may 
terminate after only a few iterations. This situation is presented in Figure [5^. 
To correct the error generated by the false bifurcation detection, child 1 is 
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Figure 4: a) Process diagram for locating child vessels at bifurcation point, b) Detection 
of the second child vessel by searching for the optimally fitting model cylinder. The search 
space is limited by the orientations of the parent vessel and child 1. 



erased and child 2 continues to track as a continuation of the parent ves- 
sel. AdditionaUy, a leak detector was implemented to remove leaks into the 
nearby mediastinum (heart and great vessels), illustrated in Figure [5}3. A 
leak is detected when a model cylinder is fit with a radius over 150% of the 
radius of the model cylinder beginning at one full cylinder length upstream. 

2.5. Sparse Surface Metric 

We developed a novel surface-based method called the sparse surface 
(SS) metric for validating the results of our pulmonary artery segmenta- 
tion method. The goals in developing a new validation criteria to be applied 
to pulmonary artery segmentation were three-fold: 

1. To quantify the number of vessels segmented (i.e. extent of segmenta- 
tion). 

2. To use the same markings to evaluate the error per vessel (i.e. the 
accuracy of segmentation). 

3. To minimize the number of manual ground truth markings (points 
defining true object geometry) required from the human user. 
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Figure 5: False positive bifurcation detection and subsequent error correcting, a) The 
child vessel does not find any strong matches, indicating that the bifurcation detector 
returned a false positive, b) The vessel has leaked into the mediastium, but is removed as 
the optimized radius grows too quickly. 



The SS method accomplishes the three listed goals by defining a mapping 
from ground truth points marked at the surface of vessels to the surface of 
the segmentation. By mapping from the ground truth to the segmentation, 
ground truth can be marked only sparsely - there is no penalty incurred 
by segmentation surface points that are not close to ground truth points. 
The full geometry of the vascular tree can be captured by marking a limited 
number of surface points for each vessel, allowing for the calculation of the 
extent of the segmentation. The following section formally defines the SS 
metric. 

2.5.1. Formal definition 

Let V be the set of voxels of the segmented volume. The boundary of V 
is the set dV of voxels that are 6-connected to a background voxel. The set 
S of surface points of a segmented volume V is then defined as: 

S = {x\x G dV} = {si, 52, Sk} (5) 

Let G be the set of ground-truth points, where the members gi were indi- 
vidually identified by manual markings. The set G contains n subsets On 
representing the ground-truth markings for n objects. The SS metric is con- 
structed from a mapping / of On to such that for every element Oi of 

On. 

f{oi) = mi = ai mm {\\{oi,Sj)\\E} (6) 

7 = 1, ...,fc 
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where the parameter describes whether the point is interior or exterior to 
the segmentation: 



. -1 if rrii ^ V (interior) 
1 otherwise (exterior) 

and II (a, 6)||^ is the positive semi-definite 3D Euchdean distance between two 
points: 

||(a,6)||^ = {a^ - h^f + (s - hyf + (a^ - h^f . (8) 

The set now holds the values of the Euclidean distance from each 
ground truth point for object n to the nearest surface point of the entire 
segmentation, with positive values corresponding to ground truth points that 
lie outside the segmentation and negative values for interior points. 

The distribution of distance elements in each of the sets is expected 
to take on a normal unimodal distribution. The ideal algorithm will return 
a narrow distribution centered around zero. For an over-segmentation, the 
distribution will be shifted towards negative distances, and for an under- 
segmentation, will be shifted towards positive distances. 

Quantitative information about the segmentation can be obtained by 
computing statistics on M^. The average signed distance (ASD) of the 
set Mn of order can be computed as: 

N 



N 

The average signed distance is a measure of the overall bias of the segmen- 
tation - a positive ASD indicates under-segmentation, and a negative ASD 
indicates over-segmentation. To evaluate the magnitude of the average error 
per truth point, the root mean square distance (RMSD) can be calculated 
for each object distance error set M^. 



nRMS 



N 9 

E^- (10) 



Additionally, the maximum amount of over- and under-segmentation can be 
found by calculating the maximum negative distance and maximum positive 
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Table 1: Statistical measurements used in the SS evaluation metric 
Statistic Measure of 



Average signed distance (ASD) 
Root mean square distance (RMSD) 
Max negative distance 
Max positive distance 



Bias of segmentation (Accuracy) 
Magnitude of error (Precision) 
Point of greatest over-segmentation 
Point of greatest under-segmentation 



distance, respectively. Table [T] gives a summary of the useful statistics in the 
SS metric. 

Given a set of surface markings for each object, the distance error statis- 
tics obtained using the SS metric can be used to achieve our goals of evalu- 
ating the number of objects segmented, and the error per object. An object 
may be classified as having been segmented if the RMS distance error is 
small, signifying that the segmentation did not stray far from the object. 
The bias, magnitude, and extreme values of the error can then be calculated 
for each segmented object. 

2.5.2. Application of the SS metric to pulmonary artery segmentation 

Ground truth points were marked on axial slices of CT scans windowed to 
soft tissue. To automatically detect segmentation errors, veins were marked 
in addition to arteries. This allowed for quantification of the main source of 
in our segmentation method. An approximately homogenous mat of points 
was created by marking vessels at roughly 4 mm intervals around the circum- 
ference of arteries for roughly every third frame, corresponding to about a 4 
mm interval in the axial direction. Each vessels was marked separately and 
given a unique label. The ASD, RMS, and maximum and minimum of the 
distances from the truth points to the nearest surface point were computed 
separately for each vessel. 

Each marked vessel was classified as either having been segmented or 
missed by measuring the RMS error per vessel. A vessel was classified as 
having been segmented if the RMS the distance was below 2.0 mm, which 
corresponds to an average error of roughly 3-4 voxels. The threshold of 2.0 
mm is an ad hoc value that represents the upper error bound for vessels visu- 
ally confirmed to be segmented. Using this method, each vessel was classified 
into one of four categories, given Table [2| for evaluating the performance of 
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the algorithm. 



Table 2: Vessel classification for evaluating algorithm performance 
Category Type of vessel Criterion 



True Positive (TP) 


Artery 


RMS 


error 


less than 2.0 mm 


False Negative (FN) 


Artery 


RMS 


error 


greater than 2.0 mm 


True Negative (TN) 


Vein 


RMS 


error 


greater than 2.0 mm 


False Positive (FP) 


Vein 


RMS 


error 


less than 2.0 mm 



2.6. Dataset 

The dataset used for this experiment consisted of 10 low-dose thoracic 
CT scans with correctly segmented and labeled airways. The scans were 
generated with an 8 slice helical CT scanner operating under a low-dose 
protocol with beam current of 40 mA and 120 kVP operating voltage. The 
slice thickness was 1.25 mm and the in-plane resolution ranged from 0.55 
mm to 0.82 mm. For each case, the left and right basal artery and all 
child branches with radii greater than 1 mm were documented with sparse 
markings. In addition, all veins with radii greater than 1 mm surrounding 
the basal arterial tree were documented. A total of 210 arteries and 205 veins 
were marked across the ten cases. 

3. Results 

3.1. Training and optimization 

Three documented cases were selected for automated algorithm training 
and parameter optimization. Table [3] summarizes the parameters that were 
optimized. 



Table 3: 


Parameter sets used in alg 


;orithm training 


Parameter 


Range Step Size 


Optimized value 


ho 


5-30 mm 5 mm 


15 mm 


^step 


0.10 - 0.30 0.05 


0.20 


^radius 


0.75 - 1.00 0.05 


0.90 
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The sparse surface metric was used to automatically calculate a training 
score S for maximization: 

S = TP + TN-FP-FN (11) 



Optimization ^radius Optimization A^^^p Optimization 

— ♦— Specificity --(■■-• Sensitivity --a-- Training Score — ♦ — Specificity • Sensitivity Training Score — ♦— Specificity Sensitivity —a— Training Score 




5 10 15 20 25 30 0.75 0.8 0.85 0.9 0.95 1 0.1 0.15 0.2 0.25 0.3 



Cylinder height (mm) Radius change threshold Incremental step size 

Figure 6: Results of the vessel tracking method for different sets of parameters. The 
optimal parameters were found by maximizing the training score. The optimal value for 
the radius change threshold (^radius reflects the approximate physiological value. 

The parameter that had the largest effect on the outcome was the cylin- 
der height ho. Visual inspection of the training results indicate that short 
cylinders are highly sensitive to noise and are unable to extract the vessel 
structure. Conversely, cylinders that were too large were not sensitive enough 
to detect most bifurcations. Similarly, incremental step sizes above the opti- 
mal value of O.2O/10 were too coarse to detect many bifurcations, decreasing 
the sensitivity. 

The radius change threshold parameter ^radius was optimized to a value 
of 0.90, which gave a marked improvement in algorithm performance over 
a much weaker value of 1.0. Limiting the search for child vessels to geome- 
tries where the radius abruptly changed improved overall vessel identification 
performance. 

3.2. Vessel tracking 

Of the 20 arterial trees to be segmented across the ten cases (10 left and 
10 right arterial trees), 19 were successfully seeded as the arteries entered 
the lungs. The vessel trackers stemming from the automated seed point 
segmented a total of 134 of the 210 manually identified arteries to within 
an RMS distance error of less than 2.0 mm, yielding an overall sensitivity 
of 64%. 20 of 205 veins were erroneously segmented as arteries, yielding 
an overall specificity for artery segmentation of 90%. Figure [7| shows the 
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Figure 7: a) Axial slice from a section of the inferior lung. Vessels in red indicate the seg- 
mentation result, b) Three dimensional rendering of the pulmonary arterial tree generated 
by our automated method. 

segmentation results, both overlaid on an axial CT slice and as a 3D surface 
representation. 

In most cases, the tracker was able to separate the arteries from the veins 
by searching for cylindrical objects in a limited size and orientation search 
space defined by the bifurcation model. While this approach was able to 
distinguish nearly identical vessels based only on their relationship to other 
arteries, it is unable to locate vessels that fall outside of the search space. 
Thus, arteries were missed if they originated at apparent trifurcations, or 
had unusually small radii relative to their sister branches. Relaxing the 
bifurcation model results in few missed segments, but more false positives 
(i.e. segmented veins). 

Figure [8] shows the results of the weakest and strongest segmentations, 
overlaid with truth markings for the arteries. In the case of the strong seg- 
mentation, the algorithm was able to correctly seed the pulmonary artery, 
and separate and segment the arterial tree to four generations. Additionally, 
the true vessel surface matched the segmented surface very closely. In the 
case of the weak segmentation, a segmental child branch failed to be detected, 
which resulted in several missed arteries. The vessel tracker also leaked onto 
the adjacent venous tree, returning false-positives (non- arteries) . 

Using the SS metric, we computed an overall under-segmentation bias 
of 0.15 mm, roughly 1/5 of a voxel. The RMS distance error between the 
manual markings and the segmentation surface was 0.63 mm, less than one 
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voxel. Among identified vessels, the greatest under-segmentation error was 
5.25 mm, which corresponded to a vessel segmentation that terminated before 
the entire true vessel was segmented. The greatest over-segmentation was 
3.14 mm. This point was found near the lung hilum, where the arteries 
and veins are closely packed, resulting in minor over-segmentation. Figure |9] 
shows the overall distribution of SS metric distance errors for the segmented 
vessels. 



Poor segmentation Strong segmentation 




Figure 8: 3D representation of the artery segmentation results, overlaid on the manually 
placed SS ground truth points, shown for a strong segmentation and an unusually weak 
segmentation. In the weak case, failure of the bifurcation detector early in the tracking 
resulted in missing several true arteries (false- negatives). 



3.3. Sparse surface metric 

To validate the use of sparse surface markings, we marked a single case 
with both dense and sparse coverage, and applied the SS metric to each truth 
set for validating the segmentation. The dense manual markings, shown in 
Figure [TO^ overlaid on a small section of the segmentation result, represent 
the most ideal coverage of the true vessel surface. In total, 4780 points were 
manually marked to completely cover the pulmonary arteries and veins of 
a single case. Figure [T0| 3 shows the same vessel covered with sparse truth 
markings. For sparse coverage, 302 points were applied to cover the same 
arteries and veins marked under dense coverage. We then asked the ques- 
tion: Under the SS metric, does sparse coverage yield the same validation 
information as dense coverage? 
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Figure 9: Distribution of SS metric distance errors for correctly identified vessels using 
sparse markings across 10 cases. The values represent the Euclidean distance between the 
manual truth points and the nearest surface point. The averaged signed distance error 
was -0.15 mm (1/5 voxel). 



The results of the SS metric distance errors are shown in Figure [TOl The 
ASD error calculated using the dense coverage was 0.29 mm, while the ASD 
error using the sparse coverage was 0.27 mm. Using an unpaired t-test, the 
95% confidence interval for the diflFerence in means between the two samples 
was (-0.07, 0.10). Thus, we find that under the SS metric, sparse marking of 
the ground truth is equivalent to more labor-intensive dense markings. 



4. Discussion 

Here we have presented a method for 1) automatically locating seed points 
within pulmonary arteries, 2) tracking the seeded arteries into the lungs 
while maintaing the separation between the arterial and venous trees, and 
3) quantitatively validating the results using a novel evaluation metric. The 
vessel tracking begins near the hilum of the left and right lungs, and uses 
a model-based bifurcation detector to segment the child vessels. In using 
a model-based method, we were able to maintain separation between the 
isointensity arterial and venous trees (90% specificity for arteries) based on 
predictable branching geometries. The segmentation was performed in binary 
space, which reduced both the signal from the vessels and the noise inherent 
to low dose scans. The arteries were separated from the veins through vessel 
morphology and relative anatomical location alone. This approach is best 
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Figure 10: Manual SS markings overlaid on the left main pulmonary artery and its two 
child branches 

suited for scans without contrast enhancement, where vessel edge information 
is lacking. 

The vessel tracking algorithm was able to successfully segment vessels 
of greatly different radii. By fitting the size of the model cylinders to the 
vessel tree, the tracker was able to segment both the large parent arteries 
of the inferior lobes and the small subsegmental arteries deep in the lung 
parenchyma. However, as the vessel radii approached the smallest possible 
cylinder, the tracker was unable to detect bifurcations. The method was 
strongest in the larger vessels where Murray's law of vessel bifurcation is 
more closely followed. As vessel sizes decrease, blood fiow becomes more 
non-Newtonian, and child vessel sizes become larger relative to the parent. 
Additionally, the smaller vessels deep in the lung parenchyma were more 
prone to obfuscation by noise because of their small size. 
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In most cases, the tracker was able to separate the arteries from the veins 
by searching for cyhndrical objects in a limited size and orientation search 
space defined by the bifurcation model. While this approach was able to 
distinguish nearly identical vessels based only on their relationship to other 
arteries, it is unable to locate vessels that fall outside of the search space. 
Thus, arteries were missed if they originated at apparent trifurcations, or 
had unusually small radii relative to their sister branches. Relaxing the 
bifurcation model results in few missed segments, but more false positives 
(i.e. segmented veins). 

The main source of error in vessel segmentation occurred at the left and 
right inferior pulmonary veins. In cases where the arteries were highly con- 
flated with the veins and airway walls, due either to high noise or unusual 
anatomy, the vessel tracker leaked onto the pulmonary veins and tracked 
them into the parenchyma. One possible method for reducing the false posi- 
tive rate would be to initiate a vessel tracker at the main inferior pulmonary 
arteries, and exclude from the arterial segmentation vessels that were highly 
likely to belong to the venous tree. 

Automated detection of parent vessels in the left and right upper lobes 
is the next appropriate step to take in the future. In this work, only the 
largest arteries to enter the lung parenchyma were identified. Such a method 
could employ a similar ROI construction or use another algorithm to locate 
the upper lobe arteries. The sensitivity of the algorithm could be improved 
by developing a trifurcation detector, or developing a method for identifying 
small subsegmental arteries missed by the tracker. Additionally, the method 
could be paired with a global Hessian matrix filter to amplify the signal from 
vessels and reduce tracking error. 

Validation of complex structures is a challenging task, and a novel surface- 
based ground truth marking system and validation metric were defined for use 
with complex structures. Sparse surface markings capture significant vessel 
detail while require fewer markings that traditional dense surface markings. 
This allows for validation on a greater number of cases given the same amount 
of effort. The sparse surface markings were shown to be an accurate sampling 
of the true surface, justifying the use of spares markings over dense markings. 

5. Conclusion 

Pulmonary artery segmentation has been suggested as a method to im- 
prove diagnosis of both pulmonary embolism and pulmonary arterial hyper- 
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tension. Although much work has focused on segmenting the arteries and 
veins as one structure, very httle previous work has been performed to distin- 
guish the two trees. In this work, a method was developed for automatically 
segmenting the pulmonary arteries in low-dose CT, and was applied to the 
basal pulmonary arteries and their child vessels. 
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